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Abstract 

We study Ua(1) anomaly of non-perturbative QCD in the maximally 
abelian gauge on SU(2) lattice. The existence of the strong correlation be¬ 
tween QCD-monopoles and instantons in the abelian gauge is shown by both 
analytic and numerical works including lattice simulations. These results 
bring us a conjecture that the Ua(1) symmetry would be explicitly broken in 
the background fields dominated by QCD-monopoles. We find an evidence for 
our conjecture by measuring the chiral-asymmetric zero modes of the Dirac 
operator in the backgrounds of QCD-monopoles. 
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As for the appearance of color-magnetic monopoles in SU(Ac) gauge theory, ’t Hooft 
proposed an interesting idea of the abelian gauge hxing [Q], which is dehned by the gauge 
transformation in the coset space of the gauge group in order to £x the gauge degrees of 
freedom up to the maximally abelian subgroup. In the abelian gauge, point-like singularities 
in the three-dimensional space under the maximally abelian subgroup can be identihed 
as color-magnetic monopoles |]l|, which will be called QCD-monopoles hereafter. In other 
words, QCD-monopoIes originate from the same topological nature as the’t Hooft-Polyakov 
monopoles, which correspond to the homotopy group 7r2(SU(A"c)/U(l)'^''“^) = p. 

However, QCD-monopoles have no condition to exist as classically stable and/or hnite- 
energy solutions, unlike the ’t Hooft-Polyakov monopoles are the explicit solutions of the 
held equation p|. One thus conjectures that the realistic model of QCD vacuum would 
be characterized by the highly quantum feature of QCD-monopoles, e.g. its condensation, 
rather than the classical one @]. Then the dual Meissner effect, which yields the exclusion 
of the color-electric helds, must be realized 0. The recent lattice QCD simulations ^- 
0 support this conjecture that QCD-monopoles play a crucial role on color conhnement 
through their condensation^]. 

As well known, QCD has also classical and non-trivial gauge conhgurations, i.e. instan- 
tons as topological defects in the Euclidean space corresponding to the homotopy group 
7r3(SU(Ac)) = Zoo It seems that instantons and QCD-monopoles are thought to be hardly 
related to each other since these topological objects appear from different non-trivial homo¬ 
topy groups. However, the recent analytical works have demonstrated the QCD-monopole 
as a classically stable solution in the background helds of the instanton conhguration using 
the abelian gauge hxing 


-2|. Furthermore, the several lattice QCD simulations have 


shown the existence of the strong correlation between instantons and QCD-monopoles in 


iQCD-monopole condensation is characterized by the presence of the long and tangled monopole 
trajectories in the four-dimensional space R^ and can be interpreted as the Kosterlitz-Thouless 
type phase transition p. 
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the highly quantum vacuum [p,|T^,|T^,|^ as well as the semi-classical vacuum 

Here, we remind that instantons are important topological objects in QCD relating to 
the Ua(1) problem m- It is well known that the index of the massless Dirac operator in 
the instanton background helds is equal to the Pontryagin index, i.e. the topological charge 

Q: 


n+ — = Q 


( 1 ) 


where n+ (n_) is the number of zero modes with the positive (negative) chirality. The 
previous relation is well known as the Atiyah-Singer index theorem. The index of the 
Dirac operator corresponds to the number of chiral-asymmetric zero modes that yield the 
chiral anomaly in the global Ua(1) symmetry. By the existence of this anomaly, the Ua(1) 
symmetry is regarded as explicitly broken at the quantum level. This mechanism plays an 


essential role on the resolution of the Ua(1) problem 


In this paper, we aim to reexamine the relation between QCD-monopoles and instantons 
through Ua(1) anomaly. By using the Monte Carlo simulation on SU(2) lattice, we measure 
the topological charge and the zero eigenvalues of the Dirac operator in both the “monopole 
dominating” and “monopole absent” background helds. 

The Maximally Abelian (MA) gauge hxing Jl| was advocated by ’t Hooft in order to 
dehne the magnetic monopole in the renormalizable and the Lorentz invariant way in the 
continuum, ± = 0 where ± In the lattice formulation this 

gauge hxing is expressed by diagonalizing the following operator X{n) through the gauge 
transformation U^{n) V{n)U^{n)V\n + fl), 


where U^{n) are link variables. While X{n) = X^{n)ai -|- X^(77,)cr2 -|- X‘^{n)a^, this gauge 
hxing means that the oh-diagonal elements are locally minimized in every sites by the 
gauge transformation, i.e. X^{n) = X‘^{n) = 0. Instead of this procedure, the gauge 
transformation is actually carried out by maximizing the gauge dependent variable R 0, 
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( 3 ) 


n, /i 

Maximizing R is equivalent to making X (s) diagonal at all sites. 

Once the gauge transformation is done by the above procedure, we factorize the SU(2) 
link variable U^{n) into the abelian link variable u^{n) and off-diagonal part Mn{n) ^ as 

Uf,{n) = M^{n) ■ Uf,{n) , (4) 

where 


M^(n) = exp{i(T 36 '^(n)} , 

(5) 

M^(n) = exp{i(TiC'^(n) i(T20^(n)} . 

(6) 


Here, 0^{n) is the U(l) gauge held and C'^(n) and C‘^{n) correspond to charged matter helds 
under a residual U(l) gauge transformation. Performing the U(l) gauge transformation on 
the original SU(2) link variable, u^{n) and M^{n) are transformed IQ as 

Uf,{n) —> %{n) = d{n)Uf,{n)d\n + fi) , (7) 

M^{n) — ^ M'^{n) = d{n)M^{n)d\n) , ( 8 ) 

where d{n) = exp{2(j393(n)}. In this way, u^{n) and M^{n) behave like an abelian ‘photon’ 
and an adjoint ‘matter’ held, respectively. 

Our next task is to look for the magnetic monopole in terms of the U(l) variables. We 
consider the product of U(l) link variables around an elementary plaquette, 

u^,y{n) = Uf,{n)u^{n R fi)u\{n -t- i>)ul{n) = , (9) 

where the abelian held strength 6^y{n) = 6y{n + jl) — 6y{n) — 6 ^{n + v) + 0^{n). It should 
be noted that the U(l) plaquette variable is a multiple valued function as the abelian held 
strength due to the compactness of the residual U(l) gauge group. Then we can divide the 
abelian held strength into two parts as 

dfiu = dfiv + 27riV^i, , (10) 
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where is the regular part dehned in —tt < 9^^, < vr and N^i, G Z is the modulo 27r 
of Here, it is known that the SU(2) link variable behaves as ~ in the MA 
gauge p. In this sense, = \s^ypaNp„ corresponds to the Dirac string following the 

DeGrand-Toussaint’s dehnition in the compact QED |]^. Then, monopole currents kp{n) 
are identihed as topological conserved currents dehned by kp{n) = dy*Npy{n + jl) [Q. 

The next aim is to extract the contribution of the monopole dominated part from the 
abelian link variable. First, we dehne the two abelian gauge helds using two parts of the 
abelian held strength [T^ as below. 


H “ rn)dx9\p{m) , 

m 

= 27r ^ G{n - m)dxNxp{m) , 


( 11 ) 

( 12 ) 


where G{n — m) is the lattice Coulomb propagator. 6^, which is called ‘regular part’, is 
composed of the regular part of the abelian held strength 9py H- On the other hand, 
which is called ‘singular part’, is composed of the Dirac string part of the the abelian held 
strength . It is noted that the sum of two values is the original abelian gauge held in the 


Landau gauge, dp9p{n) = 0 


+ C(") = Y,0{n- m)SA»Ap(m) = , 


(13) 


where a superscript L denotes the Landau gauge. This procedure could correspond to 
dividing a ‘regular photon’ part from a ‘monopole’ part . 


The corresponding SU(2) variables are reconstructed from 6^^ and 6^^ by multiplying 
the oh-diagonal factor Mp | 
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U^p\n) ^ Mp{n) , (14) 

^ Mp{n) exp{^a30°^(n)} . (15) 

Here, the fact that Up = Up^ ■ = Up^ ■ in the Landau gauge should be kept in mind. 

Thus, we shall use as the ‘monopole dominating’ SU(2) link variable and as the 
‘monopole absent’ SU(2) link variable [§- [|l0l . 
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Next, we see how instantons are defined in the lattice formulation of QCD. Of course, 
in discretised space-time, we lose inherently the topology in the strict mathematical sense. 
Nevertheless, we expect that the topological character could be neatly identified near the 
continuum limit, since the variation of fields becomes smoother than the size of the lattice 
spacing [^. We use the simplest expression for the topological charge 

1 


as 


Ql — 


327r2 


{t/^j^ (u) f/pg- (u) j- , 


(16) 


where U^y{n) is the plaquette variable. In the naive continuum limit, Tt{U^ ij{n)Upain)} is 
reduced to "lT{a^g^GpyGpa + 0{a^)} |^. The value Ql has not only O(a^) corrections, but 
also renormalized multiplicative corrections of 0(a°) Consequently, Ql is not an integer 
except for the continuum limit. However, it is known that the topological feature of Ql can 
be extracted by removing the short wavelength fluctuations of the field configuration by 


using the following procedure. We use the cooling method |^, in which each link variable 
Up is replaced by 


Up{n) = c'^U^{n)Up{n +i))Ul{n +fi) . 


(17) 


where a factor c ensures that a new link variable Up is the element of SU(2) group. After 
this procedure, which is called a cooling sweep, the action is reduced; S[Up\ < S[Up\. As a 
consequence, the field configuration locally smoothened. In Fig. [I(a)f]l(c)| , the cooling curves 
for the topological charge; Ql, the integral of the absolute value of the topological density; 
JqQ and the action divided by Svr; S are shown as typical examples. 

In order to examine the eigenvalue of the Dirac operator on the lattice, we adopt the 
Wilson fermion operator |^ ; 


D{n, 171] U) 


= 


[(Tw - %)U^(n)Sn+^^ni + (^w + “ A)' 5 n-/i, 


(18) 


iThis quantity is defined as Iq = 3 ^ ^Aii'po-|Tr{C/^,y(n)[/pcr(u)}| corresponding to the total 

number of topological pseudoparticles (instantons and anti-instantons). 
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where r„ is the Wilson parameter. We see that the choice = 1 is quite special since 1 ± 7 ^ 
are orthogonal projection operators. Thus, we use = 1 hereafter. In the naive argument, 
the Wilson fermion does not have the chiral symmetry due to the Wilson term. However, 
the effect of chiral symmetry breaking was systematically examined through the chiral Ward 
identities and the Wilson term is necessary to maintain the axial vector anomaly 


It is known that in the strong coupling region the ordinary mass term and the Wilson term 
cancel out in the pseudo-scalar mass, which is symbolically called the pion mass, at some 
K = In the strong coupling limit, the pion mass is zero at Kc(/ 3 —^ 0) ~ | |^. In 

the weak coupling regime, perturbative calculations indicate that the mass of the fermion 
becomes equal to zero along the line k = Hc{P), which ends at Kc{/3 —>• cxo) ~ |. In fact, 
many Monte Carlo simulations of lattice QCD with the Wilson fermion have shown the 
existence of such a critical line Kc{P) in the k — j3 plane, where the pion mass vanishes. The 
partial symmetry restoration would be realized near the vicinity of the critical line. 

The Dirac operator D dehned by the Wilson fermion loses a feature as the hermite 
operator owing to the discretization of the space-time. However, we can easily see that the 
operator 75 D or D^D is a hermite matrix. Then, we consider the eigenvalue problems for 
each operator by using the Lanczos algorithm. To discriminate quantum fluctuations among 
eigenvalues easily, we use the link variables smoothening the short wavelength with the 
cooling method in the same way as calculations of the topological charge The presence 
of chiral-asymmetric zero modes can be found by its sign change through the variation of 


the hopping parameter k [^. It is actually true that such a eigenvalue spectrum of 75 D 
inherently coincides with the operator D’s one sX k = Kc, where the eigenvalue spectrum 
crosses a zero line. 

We measure two sets of quantities, he. the eigenvalue spectrum of 75 D and the topo¬ 
logical charge Ql, by using ‘monopole dominating’ SU(2) link variable and ‘photon 
dominating (monopole absent)’ SU(2) link variable in the MA gauge on an 8 ^ lattice 
with fX = 2.4. Fig.|2(a)[-|2(h)| show the low-lying spectra of 75 D at various values of k for 
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8 configurations with 20 cooling sweeps in the quenched approximation. It is clear that 
there exist chiral-asymmetric zero modes around k ~ 0.132 in the background helds of the 
‘monopole dominating’ part. Furthermore, all 8 conhgurations in the background helds dom¬ 
inated by monopoles hold an analogue of “the Atiyah-Singer index theorem”; ?7,_|_ — n_ ~ Ql 
where n_|_(n_) is the number of zero modes with positive (negative) “chirality” dehned by 
the sign of its eigenvalue in the limit k | Kc- On the other hand, Fig. p(a)| - |3(h)| show that 
no existence of chiral-asymmetric zero modes is found in background helds of the ‘photon 
dominating (monopole absent)’ part in all 8 conhgurations. In this case, each topological 
charge is also equal to zero §]- p. By using 50 conhgurations, we also examine the eigen¬ 
value of a positive-dehnite hermitian operator associated with the Dirac operator D, 
in the ‘monopole dominating’ background helds and the ‘photon dominating’ background 
helds. There certainly exist almost zero modes in the ‘monopole dominating’ background, 
though we can not hud the corresponding zero modes in the ‘photon dominating’ back¬ 
ground. Therefore, in the ‘monopole dominating’ helds, the explicit breaking of the Ua(1) 
symmetry occurs due to the existence of the chiral-asymmetric zero modes. 

In conclusion, we have investigated the eigenvalue problems for the Dirac operator, which 
is dehned by the Wilson fermion, in the background helds of the ‘monopole dominating’ 
(Ds) part and the ‘photon dominating (monopole absent)’ (Ph) part by using the SU(2) 
lattice with 8^ and (3 = 2.4. In only the background helds dominated by QCD-monopoles, 
the explicit breaking of the Ua(1) symmetry occurs due to the existence of the chiral- 
asymmetric zero modes. We have found the monopole dominance for the Ua(1) anomaly, 
and also conhrmed that “the Atiyah-Singer index theorem” is satished in the backgrounds 
of QCD-monopoles. 
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FIGURE CAPTIONS 


Fig.l The cooling curves for Ql, Iq and S are examined in the ’monopole dom¬ 
inating’(Ds) part and the ’photon dominating (monopole absent)’(Ph) 
part on an 8 ^ lattice with (3 = 2.4. We show typical examples in the case 
of (a) (5l(Ds) 7 ^ 0, (b) (^^(Ds) = 0 and (c) = 0. /q(Ds) tends 

to remain hnite during the cooling process. On the other hand, /Q(Ph) 
quickly vanish by only less than 5 cooling sweeps. Therefore, topological 
pseudoparticles seem unable to live in the Ph part, but only survive in 
the Ds part in the abelian gauge. 

Fig.2 The eigenvalue spectrum of 75 D in the ‘monopole dominating’ back¬ 
ground helds (Ds part) as a function of the hopping parameter k at 20 
cooling sweeps on an 8 ^ lattice with (3 = 2.4. The chiral-asymmetric zero 
modes are found in several conhgurations, where the topological charges 
have nonzero value. 

Fig.3 The eigenvalue spectrum of 75 D in the ‘photon dominating (monopole 
absent)’ background fields (Ph part) as a function of the hopping pa¬ 
rameter K at 20 cooling sweeps on an 8 ^ lattice with (3 = 2.4. Both the 
chiral-asymmetric zero modes and the nonzero topological charges are 
not found in each conhguration. 
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